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Abstract 

In controlling distributed parameter systems it is often desirable to obtain low-order, finite- 
dimensional controllers in order to minimize real-time computational requirements. Standard ap- 
proaches to this problem employ model/controller reduction techniques in conjunction with LQG 
theory. In this paper we consider the finite-dimensional approximation of the infinite-dimensional 
Bernstein/Hyland optimal projection theory. Our approach yields fixed-finite-order controllers 
which are optimal with respect to high-order, approximating, finite-dimensional plant models. We 
illustrate the technique by computing a sequence of first-order controllers for one-dimensional, 
single-input/single-output, parabolic (heat/diffusion) and hereditary systems using spline-based, 
Ritz- Galerkin, finite element approximation. Our numerical studies indicate convergence of the 
feedback gains with less than 2% performance degradation over full-order LQG controllers for the 
parabolic system and 10% degradation for the hereditary system. 
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1. Introduction 


Approximation methods for the optimal control of distributed parameter systems have been 
widely studied. In particular, the approach taken in [1-12] involves approximating the original 
distributed parameter system by a sequence of finite-dimensional systems and then using finite- 
dimensional control-design techniques to obtain a sequence of approximating, sub-optimal control 
laws, observers, or compensators. Furthermore, in these treatments it was demonstrated that if 
the open-loop system is approximated appropriately, then it is possible to guarantee convergence 
of the sequence of sub-optimal controllers, observers, or compensators, respectively, to the optimal 
controller, observer, or compensator for the original infinite-dimensional system. In addition, it 
can also be shown that when the approximating sub-optimal control laws or estimators are applied 
to the original system, near-optimal performance can frequently be obtained. These ideas were 
pursued in the context of both open- and closed-loop control, in both continuous and discrete- time, 
and for both full-state-feedback control and LQG (i.e., Kalman-filter-based) state estimation and 
compensation. 

In practical situations, however, it is often of interest to obtain the simplest (i.e., the lowest 
order) controller which provides a given, desired feedback performance. This is usually achieved 
in one of two ways. Either the plant approximation order is reduced prior to controller design, or, 
alternatively, reduction techniques are applied to a given high-order control law. Unfortunately, 
the former approach may result in undesirable spillover effects while the latter may yield low-order 
controllers of low authority which perform unacceptably. In fact, with the second approach, this 
may occur even when a suitable controller is known to exist. For example, as is shown in [13], 
controller reduction techniques may even destabilize the closed-loop system. 

A third, more direct approach involves fixing the controller order a priori, and then optimizing 
a performance criterion over the class of fixed-order controllers. In a finite-dimensional setting, a 
set of necessary conditions in the form of four coupled matrix equations (as a direct extension of 
the pair of the separated Riccati equations of LQG theory) which characterize the optimal fixed- 
order compensator was derived in [14]. These four equations are coupled via an oblique projection 
(idempotent) matrix. In the full-order case, this projection becomes the identity thus effectively 
eliminating the additional two equations, and the necessary conditions reduce to the standard LQG 
Riccati equations. 
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The notion that this direct (i.e., fixed-finite-order) approach can be applied to distributed 
parameter systems was first suggested by Johnson in [15] and further developed in [16] and [17]. 
To realize such an approach, however, would require a suitable generalization of the optimality 
conditions for the finite-dimensional fixed-order theory. This result was subsequently obtained in 
[18] where the matrix optimal projection equations obtained in [14] for finite-dimensional systems 
were extended to a set of four coupled operator Riccati and Lyapunov equations characterizing 
optimal fixed-finite-order controllers for infinite-dimensional systems. 

In developing numerical schemes to actually compute fixed- finite- order compensators for 
infinite-dimensional systems, one might consider an approach wherein LQG reduction procedures 
are applied to a sequence of controllers obtained by using finite-dimensional full-order design tech- 
niques in conjunction with high-order finite-dimensional plant approximations. However, such an 
approach is unappealing for two reasons. First, since such methods are not predicated on the 
minimization of a performance index, prospects for convergence are slim. And, second, controller- 
reduction methods have not proven to be reliable in producing stabilizing compensators (see, for 
example, [13]). 

Hence, on the other hand, we develop an abstract approximation framework (and ultimately 
computational schemes) which combine the infinite-dimensional optimal projection theory of [18] 
with the approximation ideas developed in [9-12] for infinite-dimensional LQG problems. More 
precisely, our approach involves constructing a sequence of approximating finite-dimensional sub- 
spaces of the original, underlying, infinite-dimensional Hilbert state space along with correspond- 
ing sequences of bounded linear operators which approximate the given input, output, and system 
operators. Then, by choosing bases for these approximating subspaces and applying the finite- 
dimensional optimal projection theory, a sequence of matrix equations characterizing a sequence 
of approximating optimal, fixed-finite-order compensators for the distributed system is obtained. 
Finally, numerical techniques for solving the matrix optimal projection equations (for example, the 
homotopic continuation algorithm described in [19] and [20]) can be used to compute the sequence 
of approximating gains. 

Our primary aim in this paper is to describe the general approach we are proposing, to discuss 
its implementation, and to demonstrate its feasibility and practicality. We offer no convergence 
arguments here, but rather reserve them for a more theoretical paper to follow. Instead, we consider 
the application of our technique to two examples. One involves the control of a one-dimensional, 
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single-input, single-output parabolic (heat/diffusion) system while the other involves a single-input 
single-output one-dimensional hereditary control system. These relatively simple examples have 
been used throughout the distributed parameter control literature to illustrate the application of 
new theories and techniques. A detailed discussion of the application of our ideas to more complex 
control systems, for instance, the vibration control of flexible structures, will also appear elsewhere. 
We use spline- based Ritz-Galerkin finite element schemes to approximate the open-loop systems 
(one for which convergence can be demonstrated in the LQG case) and present and discuss some 
of the numerical results which we have obtained using our general approximation framework. 

We now outline the remainder of the paper. In Section 2 we briefly review the infinite- 
dimensional optimal projection theory from [18], describe the approximation framework, and derive 
the corresponding equivalent matrix equations and feedback gains which characterize the approx- 
imating fixed-finite-order compensator. In Section 3 we consider the examples, construct the ap- 
proximation schemes, and discuss our numerical findings. Section 4 contains a summary and some 
concluding remarks. 

2. Optimal Projection Theory and Finite-Dimensional Approximation 

We consider the following fixed-finite-order dynamic-compensation problem. Given the infinite- 
dimensional control system 

x(t) = Ax(t) + Bu(i) + Hiw(t) (2.1) 

with measurements 

y(t) = Cx(t) + H 2 tv(t), (2.2) 

where t E [0,oo), design a finite-dimensional, n c th-order dynamic compensator 

x c (t) = A c x c (t) + B c y(t), (2.3) 

u(£) = C c x c (t ) (2.4) 

which minimizes the steady-state performance criterion 

J(A c ,B c ,C c ) = Urn JE[{R x x{t),x{t)) + u(£) T i? 2 u(i)] . (2.5) 

For convenience we denote the infinite-dimensional plant by 77; that is, 

II±{A,B,C,R U R 2} V U V 2 }. 
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Here x (£) lies in a real, separable Hilbert space X with inner product (•, •), A : Dom (A) C X — * X is 
a closed, densely defined operator which generates a C 0 semigroup {T(t) : t > 0} of bounded linear 
operators on X,B G £( IR m ,X), and C G £(X We assume that the state and measurement 
are corrupted by a white noise signal u>(t) in a real, separable Hilbert space X (see [21] or [22]), 
that Hi G £(X , X),H 2 E £(.£, IR*), R± G £{X) is (self-adjoint) nonnegative definite, and that R 2 
is an m X m (symmetric) positive-definite matrix. We define Vi = and V 2 = H 2 H^ where 

( )* denotes adjoint, and assume for convenience that HiH £ = 0 and that V 2 is positive definite. 
The compensator is assumed to be of fixed, finite order n c (i.e., x c (i) E IR n °) and that A C ,H C , 
and C c are matrices of appropriate dimension. For further details and discussion on the problem 
statement and the above assumptions, see [18]. 

We summarize here the primary result from [18] characterizing optimal fixed- finite-order con- 
trollers. For convenience define H = and H = C*Vj~ 1 C. Also let I ne and lx denote 

respectively the n c x n c identity matrix and the identity operator on X . 

Theorem 2.1. Let n c be given and suppose there exists a controllable and observable n c th- 
order dynamic compensator (A c , B cy C c ) which minimizes J given by (2.5) and for which the closed- 
loop semigroup generated by 

A F A BC c 
[B c C A c 

is uniformly exponentially stable. Then there exist nonnegative-definite operators Q,P, Q,P on X 
such that A C ,B C , and C c are given by 


A c = r(A-Q2 - ZP)G*, 

(2.6) 

B c = rQC'Vf 1 , 

(2.7) 

Cc = -Ri l B‘PG\ 

t. 

(2.8) 

where Q : Dom (A*) — ► Dom(A),P : Dom (A) — * Dom(A*),<5 : 

: X -* Dom (A), P : X -* 


Dom (A*), and G,T G C(X , IR nc ), and such that the following conditions are satisfied: 

rank Q = rank P — rank QP = n c , (2.9) 


QP = G*MF, rG* = I ne , 


( 2 . 10 ) 


for someMer^"”, 
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0 = AQ + QA* + Vi - QSQ + t x QSQt * x , (2.11) 

0 = A*P + PA + R ± - PSP + t^PSPtx, (2.12) 

0 = (A - SP)Q + Q(A - SPY + QSQ - t x QSQt * x , (2.13) 

0 = (A - Q£)*P + P(A - Q£) + PSP - t x PSPt x , (2.14) 

where 

t = G*r, r x = /r-r. 

It is shown in [18] that the factorization (2.10) for the nonnegative-definite operators Q and P 
satisfying rank QP = n c always exists and is unique except for a change of basis in It is also 
shown in [18] that 7 1 * : JR n ° — ► Dom(A*) so that the expression (2.6) is well defined. 

Equations (2.11)-(2.14) are, in general, infinite-dimensional operator equations. To actually use 
them to compute the optimal fixed-order compensator, a finite-dimensional plant approximation is 
required. For each IV = 1,2, . . . , let X N denote a finite-dimensional subspace of X and let P N : 
X — > X N denote the corresponding orthogonal projection of X onto X N . Let A N G £(X N ),B N G 
£(IR m , X N ),C N G £(X N ,m l ),R? G £{X N ), and V X N G £{X N ). We consider the system (2.6)- 
(2.14) with the plant 77 replaced by the plant 77 N given by 

n N ^{A N i B N i C N 9 R^ } R 2} V 1 N i V 2 }. 

Typically, the operators B N ,C N , R ± and Vf* are chosen as B N = P N B,C N = CP N y = P N Ri 
and Vf* = P N Vi with the requirement that P N converge strongly to the identity I x as N — ► oo. 
The operator A N is chosen so that it and its adjoint satisfy the hypotheses of the Trotter-Kato 
semigroup approximation theorem (i.e., stability and consistency, see, for example, [23]). That is, 
A n is chosen so that lim/z^oo T N {t)P N <f> = T(t)<f > , and lim/^oo T N (t)* P N <f> = T(t)*<f > , uniformly 
in t for t in bounded intervals, for each <f> G X , where T N (t) = exp(iA^), t > 0. We shall say more 
about these choices for A N 7 B N y C N > 72^, and when we address convergence questions below. 

Although with the plant 77 N equations (2.11)-(2.14) are finite dimensional, they are still oper- 
ator equations. It is their matrix equivalents which are used in computations. Unless orthonormal 
bases are chosen for the subspaces X N (which is typically not the case in practice) some care must 
be taken to obtain the appropriate matrix system. 
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For each N = 1,2, let {^*}y=i be a basis for X N and choose the standard bases for all 
Euclidean spaces. For a linear operator L with domain and range X N or any Euclidean space, 
let [. L ] denote its matrix representation with respect to the bases chosen above. Also, let <P N 
denote the A: -^-square Gram matrix corresponding to the basis {0*}y=ij that is, = (<j>? , <^), 
i,j = 1, 2 , . . . , k N . Noting that 

[(A*)*] = (<P*) _1 [A*] T **, [(B*)*] = [(C N )] = ($*)- 1 [C*] T , 

[(jf )*] = (^)" 1 [rx ]' V, [27*] = [B^BJ-^B*]*#*, . [5*] = (^^[C^fMC*], 

the matrix equivalents of the operator equations (2.11)-(2.14) become 

0 = [A*][Q*] + [Q*](£*) _1 [A*] T £* + [V*] - [Q*][r*][Q*] 

+ [rfl[Q "PVl - (^)- 1 [r*l V, (2.15) 

0 = (<?*)- 1 [A*] t 4>*[P 7V ] + [P N ] -‘[A*] + [R?] - [P*][27*][P*] 

+ («P*)- 1 [r*] T «P*[P*] - [27*][P*][rf ], (2.16) 

0 = ([A*] - [27*][P*])[Q*] + [Q*]($*) _1 ([A*] - [S N ][P N ]) T $ N 

+ [Q*][U*][g*] - [r*][g*][r*][g*](<P*) _1 - [r*] T **, (2.17) 

o = (#*)-» ([a*] - [q n ][s n ))V[P n ] + [p*]([a*] - [g*][r*]) 

+ [P*][27*][P*] - ($* ) “ 1 [r*] [P*] [27*] [P*] [r*] . (2.18) 

Therefore, if we define the k N X k N nonnegative-definite matrices 

Qo = [g*](#*) _1 , P 0 *±#*[P*] f 

Qo = [q*](<p*) -1 , P 0 N ±*"[P N ], 

v 0 * = KK^*)- 1 , < = **[*?], 

27* 4 [P^]f? 2 - 1 [B JV ] T , £ 0 * = [C*]nr l [C N ], 

we can solve the matrix optimal projection equations given in [14] corresponding to the matrix 
plant model 

n 0 N ±{[A%[B”],[C”],RZ,R 2 ,V Q ",V 2 }, 

to obtain the matrices <3*,P*,<3* an d P*. The approximating optimal n c th-order dynamic 
compensator {A*,.B*,C*} is then given by 

A* = P 0 *([A*] - Qo^o ~ Z 0 N Po N )WV> 
b? = Pc ?Qolc N ] T v 2 -\ 

C? = -i? 2 - 1 [B*] T P 0 *(^) T , 
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where G IR" oXfcN , and M N G nt" <Xn ° satisfy 

Q”P 0 n = G”M N r”, r 0 N (G ") T = I ne , 

[r w ] = (G”rr 0 N , [rf] = I k u - [r w ]. 

When an infinite-dimensional controller will suffice, C c = —R^B'P G C(X ,IR m ) and B c = 
QC^f 1 G £{R l , X) are the usual infinite-dimensional LQG controller and observer gains (see [9]). 
The operators P 7 Q G C(X) are the nonnegative-definite solutions to the two decoupled operator 
algebraic Riccati equations (2.11) and (2.12) with t and formally set to lx and 0, respectively. 
Since C c has range in IR m and B c has domain IR*, there exist vectors c c = (cj, . . . ,c^) T G X 
and b c = (&J, . . . ,6£) T G x£ =1 X such that 

[C c x] i = {c i c9 *), t = 1,2, . . . ,m, xGX 7 

and 

t 

B c y = &Jy = '^2y i b t c , y = (y 1( . . . , y t ) G IR 4 . 

t=l 

The vectors c c and b c are referred to as the optimal LQG controller and observer functional gains 
respectively. 

With regard to approximation for the full-order LQG problem, for each N = 1,2, ... we take 
n c = k N . Then it is not difficult to show that 

C?[P N X] = <c" *), xex, 

and 

B?y = (&") T y, y € ni 4 , 

where c G XjLiX N and b G x 4 =1 jf N are given by c % = {^ N )~ 1 <f> N and 6 P = (B?) T <f> N 

respectively with <f> N = ((f>[ f ,...,<f>% N ) G Xj^X 1 *. The vectors and b P are referred to as the 
approximating optimal LQG controller and observer functional gains. To compute them we need 
only solve two standard decoupled matrix algebraic Riccati equations for the k N x k N nonnegative- 
definite matrices Qq and P^. 

A rather complete convergence theory for LQG approximation can be found in [9]. Essentially, 
it is shown there that if the approximating subspaces X N are chosen so that the projections P N 


7 



converge strongly to the identity as N — ► oo, the operators A N ,B N ,C N ,R ^ , and Vf* are chosen 
as was described above, and the operators Q N and P N are uniformly bounded in N, then Q N and 
P N converge weakly to Q and P, respectively as N — ► oo. This in turn implies that C ^ — + C c , 
strongly, B ? — ► B c , weakly, — ► c c and b^ — + b c , weakly, and the closed-loop semigroup 
for the approximating optimal LQG compensator converges weakly to the closed-loop semigroup 
for the optimal infinite-dimensional LQG compensator, as N — ► oo. If, in addition, the operators 
S N (t) = T N [t)+B N C n and S N (t) = T N [t) — B^C N are uniformly exponentially stable, uniformly 
in N, then Q N — ► Q and P N — > P, strongly, C ? — ► C c and B% -* B C) in norm, c —>■ c c and 

b —> b c , strongly, and the closed-loop semigroups converge strongly, as N — * oo. If and 

are coercive, uniformly in N, then S N (t) and S N (t) will be uniformly exponentially stable. If it is 
also true that Ri and Vi are trace class and R^ P N — ► Pi and Vf* P N — * Vi in trace norm then Q 

and P are trace class and Q N P N Q and P N P N —> P in trace norm as N — » oo. 

Returning to the fixed-finite-order czLse, we note that in general the approximating optimal 
projection equations may not possess a unique solution. However, in [19] it is shown for the finite- 
dimensional case that it is possible to obtain an upper bound for the number of stabilizing solutions. 
Using topological degree theory, the following result was obtained in [19]. 


Theorem 2.2* Consider the equations (2.1l)-(2.14) with the infinite-dimensional plant 77 re- 
placed by the finite-dimensional plant 77 N . Let n u denote the dimension of the unstable subspace of 
A n and assume that n c > n u . Then in the class of nonnegative-definite operators Q N ,P N ,Q N ,P N 
on X N satisfying rank Q N = rank P N = rank Q N P N = n c , there exist at most 

fmm(k N ,m,l)-n u \ minfh N m .x 

I I ) — mini A: , m, £), 

\ n c - n u / 

1, otherwise, 

solutions of (2.11)-(2.14), each of which is stabilizing. If, in addition, the plant (A N , B N ,C N ) 
is stabilizable by an n c th-order controller, then there exists at least one stabilizing solution of 
(2-9)— (2.14) . 


Theorem 2.2 shows that while there may exist multiple solutions to the finite-dimensional op- 
timal projection equations, in practice this number can be quite small. For example, if n c > n u 
and the system is either single input (m = 1) or single output (£ = 1) then there exists at most 
one solution to (2.9)-(2.14) for the plant II N . The existence of at least one stabilizing solu- 
tion of course depends upon whether or not the plant is stabilizable by an n c th-order controller 
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(for relevant results, see [24]). Finally, while it may be possible to stabilize a plant with n c < n u , 
this case lies outside the scope of the analysis given in [19]. 


3. Examples and Numerical Results 


We first consider the one-dimensional, single-input/ single-output, parabolic (heat/diffusion) 
control system with Dirichlet boundary conditions given by 


dv 


d 2 v 


■g^{t, 1 l) = a g^(t,r}) + b(i])u{t) + h 1 (ii)w 1 (t,T)) ) 0 < t) < 1, t> 0, 

v(t, 0) = 0, v(i, 1) = 0, t > 0, 

y(t)= [ c{ri)v{t,ri)d.r] + h‘ l w i {t), t> 0, 

Jo 

where a > 0, and &(•) and c(*) are given by 


and 


[ 0, elsewhere, 

4 ») = { ^ I s 1,1 

[0 } elsewhere, 


(3.1) 

(3.2) 

(3.3) 


with 0 < /?i < /? 2 < 1 and 0 < 71 < 72 < 1. In (3.1) and (3.3), /i(-) G Loo(0, 1), u>i(£, •) G 
L 2(0, 1), a.a. t G [0,oo), (see [22], p. 314), is a nonzero constant and u/2(*) is unit-intensity white 
noise. 


To rewrite (3.1)-(3.3) in the form (2.1), (2.2), in the usual way we take X = £2(0, 1) endowed 
with the standard L 2 inner product, let x(t) = v(t,*),t > 0, define A : Dorn (A) C X — ► X by 
A(j> = aD 2 <f> for <f> G Dorn A = R 2 (0,l) ni?^(0,l), and define B G ^(IR 1 ,!) and C G C{X ^JR 1 ) by 
Bu = 6(-)u for u G IR 1 , and C<f> = c{r))<f>{r})dr { j, for <f> G -Z^2 (0» 1), respectively. Furthermore, let 

X = Zr2 (0, 1) X IR, set ty(i) = (tx^x (^ , •) , W2 (*)) € X-> and define Hi G £(X,X) and H 2 G £(JC',IR 1 ) 
by Hiz = hi{-)zi and H 2 z = h 2 z 2 for z = (zi,z 2 ) £ X. 


It is well known (see, for example, [23]) that A is closed, densely defined, and negative definite. 
Furthermore, A is the infinitesimal generator of a uniformly exponentially stable, analytic (abstract 
parabolic) semigroup {T(t) : t > 0} of bounded, self-adjoint linear operators on X . 

We consider linear spline-based Ritz-Galerkin approximation for the open-loop system. For 
each N = 2,3, . . ., let be the linear spline (“hat”) functions defined on the interval [0,1] 
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with respect to the uniform partition {0? jj, - . . , 1} , i.e., 


f Nrj-j + l, 


4? in) = 


ji + l- Nrj, 

U 


*> e ftvT n)> 

e [iv’ *i v~)» 
elsewhere on [0,1], 


j — 1,2, . . . y N — 1. Set X N = span and note that k N = dimX^ = N — 1, and 

X N C ffo(0, 1) for all N, If P N : X — > X^ denotes the orthogonal projection of X onto X N , 
then standard convergence estimates for interpolatory splines (see [25]) can be used to show that 
lim j tf_> 00 P N <f> = <f> in £ 2 ( 0 , 1) for (f> G £ 2 ( 0 , 1). 


There are two equivalent ways to obtain an operator representation for the usual Ritz-Galerkin 
approximation to A. First, A can be extended to a bounded linear operator from #q( 0, 1) onto its 
dual, if” 1 (0,1), via 

[A<}>){4 >) = € *£(0,1). (3.4) 

Since X N C #o(0,l) for all N = 2,3, . . . , we define A N G ^(X^) by A N <f> N = A<f> N ,<f> N G X^, 
with A^ G #“ 1 (0, 1) considered to be a linear functional on X^. From the Riesz Representation 
theorem we obtain A N <f> N = t p N where 1 p N is that element in X^ which satisfies (A iV ^ iV )(x iNr ) = 

-a(D(f> N ,Dx N ) = {^ N ,X N )- 

Alternatively and equivalently, by using the fact that A is self-adjoint, we can define A N as 
follows. Let Pf* : #q(0, l) — * X N denote the orthogonal projection of the Hilbert space # 0 ( 0 , 1 ) 
onto X*. Using the definition (3.4), it is not difficult to show that —A G £(#q (0> 1), # _1 (0, 1)) 
is coercive and, therefore, that A -1 : #“ 1 (0,1) — ► #q(0,1) exists and is bounded. We then define 
A n G £(X n ) to be the inverse of the operator given by (A^) -1 = Pf* A~* \x*r . 

Using either definition, it is easily argued that A N is well defined, self-adjoint, and is the 
infinitesimal generator of a uniformly exponentially stable (uniformly in N ) semigroup, T N (t) = 
exp (tA N ),t > 0, of bounded linear operators on X^. Also, using the approximation properties of 
splines, it is not difficult to show that lim^_^ 00 (A iV ) _1 / ? ^^ = A~ l <f>, (f> G X. Consequently, the 
hypotheses of the Trotter-Kato theorem (see [23]) axe satisfied and we have lim J N'_*. 0 o T N (t)P N <f> = 
T(t)<f> and limjv_*oo T N (t)*P N <f> = T(t)*<f> y <f> G X, uniformly in t for t in bounded intervals. A 
detailed discussion of the results just outlined can be found in [8]. 

We define B N = P N B and C N = CP N , from which it immediately follows that linrvN , _>oo B N = 
B and lim^^oo C N — C in norm and similarly for their adjoints. For the example we shall consider 
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here, we have chosen Ri = r^Ix^Ri ~ film, with ri,r2 > 0. Setting hi(rj) = ,0 < ij < 1, and 

h 2 = v£ with vi,U2 > 0, we obtain V\ = lx and V 2 = v 2 . We then take R^ = and 

= P N V i. For the LQG problem, the open-loop uniform exponential stability of both the 

infinite-dimensional system and the approximating systems is sufficient to conclude the strong 
convergence of the approximating Riccati operators to the solutions of the infinite-dimensional 
Riccati equations, the uniform norm convergence of the approximating controller and observer 
gains, and the strong convergence of the functional gains, as N — * oo. 

Since the basis elements are piecewise linear with respect to the uniform mesh 

{0, 1} on [0,1], the equivalent matrix representations for the operators defined above can 

be computed directly and in closed form. The Gram matrix = 1,2,...,N — 1 

is given by $ N = -“Tridiag{-, and if we define the generalized stiffness matrix & N by 

= —a(D<f>^,D(f>^) y t i j = 1,2,. ..,iV — 1, then }F N = aN Tridiag{l, — 2, 1}. It follows that 
[A n ] = (# W )- 1 <F JV ,[£ W ] = ($ N )~H N ,{C N ] = c N , with Sp? and 

cf = (c,«£f ) = ^3^7 = 1,2, . • • ,N- 1, and that Rg = r l ^ N and V 0 N = 

For our numerical study we set a = l,/?i = .75 — .03\/2,/?2 = 75 + .04\/2, 71 = .25 — .04\/2,72 = 

.25 + .03\/2, ri = t>i = 1, r*2 = t>2 = 10~ 4 , /ii (77) = 1, and used our technique to compute approxi- 
mating optimal LQG (i.e., n c = N — 1) and 1st order (i.e., n c = 1) compensators for various values 
of N. The open-loop stability of system (3.1)-(3.3) and the approximating systems imply that the 
finite-dimensional approximating optimal projection equations have a solution. Theorem 2.2 on the 
other hand, with n u = 0 and n c = 1 or n c = N — 1, implies that they have at most one solution. 
Consequently, the system of equations (2.11)-(2.14) with the plants II N admits a unique solution. 

The optimal projection equations (2.11)-(2.14) were solved using the homotopic continuation 
algorithm described in [19]. It is shown in [19] that the operation count for the algorithm is 
proportional to p(2n 3 + (m + £)n 2 + (m + £) 3 n 3 ) where p is the number of integration steps and n 
is the dimension of the finite-dimensional plant. This is competitive with the operation count for 
the Hamiltonian solution of the standard Riccati equations which is 0(16n 3 ) for LQG. Also, note 
that the computational burden for the solution of the optimal projection equations decreases with 
n c . 


Since m — t~ 1 in the LQG case, the optimal functional observer and feedback control gains 
b c and c c and the approximating gains 6 ^ and c^, are all simply L 2 functions with b ** and c ** 
elements in X N , We plot the functions b ^ and we obtained for various values of N respectively 
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In the fixed-order case with n c = 1, the compensator gains A c , B c , and C c are all scalars. Also, 
for a first-order controller there axe only two independent parameters, A c and B c C c . In Table 3.1 
below we give the values we obtained for A** and for various values of N. Once again, 

it is clear that the gains are converging as N increases. In addition, in Table 3.1 we provide the 
closed-loop costs J^qg ^fo f° r the LQG and first-order controllers. These closed-loop costs 
were evaluated using a 64th-order modal approximation to the infinite-dimensional system. For all 
values of N the performance of the fixed-order compensator was within 2% of the corresponding 
LQG controller. Thus, for example, the replacement of a 32nd-order approximating optimal LQG 
controller by an approximating optimal first-order controller will yield considerable implementation 
simplification with only minor performance degradation. Note that for the example we consider 
here, it is possible to compute the open-loop cost for the infinite-dimensional system in closed form. 
We have 


Jql — tr 





Uxritr 



T{t) 2 dt 


= v i r i 



e -2 n** 3 at dt _ 


o° . 

^>•1 y' J_ 

2 n 2 a n 2 

n = 1 


ViTi 

12a 


« .08333. 
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Finally, for comparison purposes, we tried applying balancing techniques to the LQG controllers 
to reduce their order. However, with n c = 1, such controllers were found to be destabilizing. Based 
upon the results in [13], this was not unexpected. 


N 

A" 

b n c n 

J N 

J LQG 

J?o 

4 

-687.6 

5470 

.06999 

.07014 

8 

-720.9 

5231 

.06870 

.06993 

12 

-730.9 

5182 

.06872 

.06991 

16 

-734.3 

5145 

.06874 

.06990 

20 

-738.0 

5127 

.06875 

.06990 

24 

-737.6 

5108 

.06876 

.06990 

28 

-739.8 

5109 

.06876 

.06990 

32 

-738.7 

5099 

.06877 

.06990 


Table 3.1 
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As a second example we consider the one-dimensional, single-input, single-output hereditary 
control system given by 

v(t) = a 0 v(t) + aiv{t — p) + b 0 u(t) + hiw(t) y t > 0, (3.5) 

y (*) = c 0 t;(t) + h 2 w(t ), t > 0, (3.6) 

where ao, ai, &o> Co, hi, h 2y p £ IR 1 with h 2 ± 0, and w is a unit- intensity white noise process. To 
rewrite (3.5), (3.6) in the form (2.1), (2.2), we take X = IR 1 X L 2 (— p y 0) endowed with the usual 
product space inner product, {(*7, <£),(£,*/;)) = *]£ + an( l let = (v(i),t; t ), * > 0, where 

for t > 0, v t £ L 2 (—p y 0) is given by tr t (0) = v (* + 0)> — p < 0 < 0. Define A: Dom(A) C JC — > X by 
A4> = (a o m + ai<p(—p), D(j>) for $ = {<f>{ 0),<f>) G Dom(A) = {{&$) G X:ip G H x {-p, 0),^(0) = 
f}, and let B G ^(IR 1 ,^) £ind C G £(X’,IR 1 ) be given by Bu = (6 0 «,0) and C(rj,<f>) = cqT}, 
respectively. Let X = IR 1 and define Hi G C(X, X) and H 2 G ^(XjIR 1 ) by H\z = ( hiz,0 ) and 
H 2 Z = h^z, for z G IR 1 . 

The operator A is densely defined and is the infinitesimal generator of a C 0 semigroup {T(t): t > 
0} of bounded linear operators on X with T{t)(r] y (j>) = (v(t;r] y <f>) y Vt(r) y <£)), t > 0, where v{'\r) y <f>) 
is the unique solution to (3.5) with b 0 = hi = 0, and initial conditions v(0) = r] y v 0 = (f>. We take 
Ri £ £(X) and R 2 £ ^(IR 1 ) to be Ri{rj y <j>) = (ri77,0) and R 2 u = r 2 u, respectively, with ri,r 2 > 0. 
The definitions of Hi and H 2 given above imply that Vi £ C(X) and V 2 G ^(IR 1 ) are given by 
= (^ l i r 7> 0) and V 2 z = h\z y for ( r] y <f > ) £ X and z £ IR 1 . 

We employ an approximation scheme recently proposed by Ito and Kappel in [26]. We briefly 
outline it here; a more detailed discussion can be found in [26]. For each N = 1,2, ... let xf G 
L 2 (— /?, 0) denote the characteristic function for the interval [— jp/N y —(j — 1 )p/N) y j = 1,2,. . . , N y 
and let X N be the ( N + l)-dimensional subspace of X defined by„ 

X N = span{(l,0),(0,x^),..-,(0,xjv)}- 

Let P N :X — ► X n denote the orthogonal projection of X onto X N . Let denote the linear B- 

spline functions defined on the interval [— p,0] with respect to the uniform mesh [— p y . . . , —p/N, 0}, 
and set X ^ — span {{^ (0), Then Xi is an (iV + l)-dimensional subspace of Dom(A) 
and it is not difficult to demonstrate that the restriction of P N to Xj* is a bijection onto X N . Using 
the fact that A restricted to X has range in X N , we define A N G £(X N ) by A N = A[P N )~ l , 
and set T N (t) = exp(A N t), t > 0. Noting that R(B) C X N , we take B N G ^(IR 1 , X N ) to be given 
by B N = B. Similarly, we take = R\ and Vf* = Vj. We set C N = C. 
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It is shown in [ 26 ] that P N {r},<f>) —* (ti,<f>), T N (t)P N (ri,<f>) — > T(t)(r],(f>), and 

T N (t)* P N (r},<f>) — > T(t)*{r],<f>) for (»?,^) € X as IV — » oo, uniformly in i, for i in bounded subsets 
of [ 0 ,oo). It then follows that limjv-*oo B N = B and limjv-^oo C N P N = C, in norm. 

For the LQG (full-order) problem, the optimal functional observer and feedback control gains 
b c and c c are of the form b e = and c c = (70,71) with /?o,7o G IR 1 , and /?i,7i € ^2(-p,0). 

The approximating gains are of the form b = ( 0 o , 0 ^) and = (7^ ,7^) with 0 o,lo G IR 1 and 
i'll G span {x^}yLr Since we are treating a one-dimensional example, if 60 7^ 0 , the theory 
in [ 26 ] implies that ~ * A) and 7^ — ► 70 in IR 1 , and 0 i —* 0 i, and 7^ — » 71 in Li(—p, 0), as 
JV-» 00. 

Once again, as in the first example, matrix representations for the operators A N , B N , C N , , 

and V-^ are not difficult to compute in closed form. Indeed, the (N + 1 ) X (N + 1 ) matrix repre- 

0 

0 * * 0 
0 \ *J 


0 aj 

0 * _AL 0 

p p 

0 K _N 

P P J 

We have the (IV + 1 ) x 1 matrix [B w ] = [i 0 0 . . . 0 ] T and the 1 X (N + 1 ) matrix [C^] = [co 0 ... 0 ] , 
while [R*] = rx[Af w ] and [V/^] = /i f ^ ] where the (IV + 1 ) X (IV + l) matrix [W w ] is given by 

0- 
0 

0. 

We set a 0 = ai = bo = c Q = ri = hi — p = 1 , r% = .1, and and computed 

approximating optimal LQG (i.e., n c = N + 1) and first-order (i.e., n c = 1) compensators for 


[M N \ = 


r 1 0 
0 

Lo 


sentation for the bijection P N : X{* — ► X\ is given 1 


[p n ) = 


1 0 
i i 0 
0 i i 0 


Then [A w ] = [Ff^JfP^] 1 , where 


[K n ] = 


a 0 0 

N_ _ff 
P P 

0 tL 


-1 ° 
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N = 8,16,24, and 32. The optimal LQG observer gains are given in Table 3.3 and Figure 3.3; the 
control gains are given in Table 3.4 and Figure 3.4. The first 23 open-loop poles of the system 
(see [27]) are given in Table 3.2. The approximating first-order compensator gains along with the 
corresponding and LQG closed-loop costs are given in Table 3.5 below. These costs were computed 
using an evaluation model obtained by setting N = 64.. Note that the performance of the first- 
order controllers is within 10% of the performance of the LQG controllers. Once again it is clear 
that convergence is achieved. 


1.278465 

-1.588317 ± 4.155305i 
-2.417631 ± 10.68603i 
-2.861502 ± 17.05611i 
-3.167754 ± 23.38558i 
-3.401945 ± 29.69798i 
-3.591627 ± 36.00146i 
-3.751047 ± 42.29965i 
-3.888543 ± 48.59442i 
-4.009422 ± 54.88686i 
-4.117267 ± 61.17761i 
-4.214618 ± 67.46710i 


Table 3.2 
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_N_ 8 16 24 32 

7^ -4.4213 -4.4229 -4.4233 -4.4234 

Table 3.4 
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Figure 3.3 


JV 8 16 24 32_ 

4.4213 4.4229 4.4233 4.4234 

Table 3.3 
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N 

A? 


Jlqg 

J?o 

8 

-4.835 

-16.057 

1.4042 

1.5221 

16 

-4.936 

-16.343 

1.403877 

1.5298 

24 

-4.959 

-16.378 

1.403856 

1.5309 

32 

-4.962 

-16.404 

1.403852 

1.5317 


Table 3.5 


4. Summary and Concluding Remarks 

We have proposed an approximation technique for computing optimal fixed-order compensators 
for distributed parameter systems. Our approach involves using the optimal projection theory for 
infinite-dimensional systems (which characterizes the optimal fixed-order compensator) developed 
in [18] in conjunction with finite-dimensional approximation of the infinite-dimensional plant. We 
demonstrated the feasibility of our approach with two examples wherein we used spline-based Ritz- 
Galerkin finite element schemes to compute approximating optimal first-order controllers for one- 
dimensional, singe-input/output, parabolic (heat/diffusion) and hereditary control systems. Our 
numerical studies indicate that convergence of the compensator gains is achieved and that using 
the first-order controller would lead to only minimal performance degradation over a standard LQG 
compensator while yielding significant implementation simplification. 

At this point one is led naturally to ask the question of whether or not a satisfactory convergence 
theory could be developed. We are working on this at present and expect that such a theory 
would conform closely in form and spirit to the convergence results for LQG approximation found 
in [9] and [10] and outlined in Section 2 above. We also intend to consider our approximation 
ideas in the context of discrete-time or sampled-data systems, and for continuous-time systems 
involving unbounded input and/or output (for example, boundary control systems), and systems 
with control or measurement delays, see [11], [12]). Finally, we intend to investigate the application 
of our approximation framework to other infinite-dimensional control systems, in particular the 
vibration control of flexible structures (i.e., second-order systems such as wave, beam, or plate 
equations). 
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